EXPLICIT CONSTRUCTIONS OF QUASI-MONTE CARLO RULES 
FOR THE NUMERICAL INTEGRATION OF HIGH DIMENSIONAL 

PERIODIC FUNCTIONS 

JOSEF DICK* 

Abstract. In this paper we give explicit constructions of point sets in the s dimensional unit cube 

p^ ' yielding quasi-Monte Carlo algorithms which achieve the optimal rate of convergence of the worst-case 

error for numerically integrating high dimensional periodic functions. In the classical measure Pa 

of the worst-case error introduced by Korobov the convergence is of 0{N~ mm(Q,d) {\Qg A^)"""^) for 

every even integer « > 1, where d is a parameter of the construction which can be chosen arbitrarily 

large and N is the number of quadrature points. This convergence rate is known to be best possible 

5-H ' up to some logAf factors. We prove the result for the deterministic and also a randomized setting. 

S—^ The construction is based on a suitable extension of digital (t, m, s)-nets over the finite field Z;,. 
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^H . 1. Introduction. Korobov [13] and independently Hlawka [llj introduced a 

r~| ' quadrature formula which is suited for numerically integrating high dimensional pe- 

riodic functions. More precisely, we want to approximate the high dimensional inte- 
gral J,„ ^,^ f{x) Ax (where / is assumed to be periodic with period 1 in each coordi- 
nate) by a quasi-Monte Carlo rule, i.e., an equal weight quadrature rule QN,s{f) = 
N~^J2n=o /(^")' ■where Xo,. . . ,X]y^i € [0,1]"* are the quadrature points. Specifi- 
cally, Korobov and Hlawka suggested using a quadrature rule of the form QN,g,sif) = 
Q^ , N~^ X)n=o /({"-S/^})' where for a vector of real numbers x = {xi, . . . , Xg) we define 

CN I {x} as the fractional part of each component of a;, i.e., {xj} = Xj— [x^J = Xj (mod 1) 

^L) • and where g G Z* is an integer vector. The quadrature rule QN,g,s is called lattice 

rule and g is called the generating vector (of the lattice rule). The monographs 



■"nJ" . [m [H [HI [27] deal partly or entirely with the approximation of such integrals. (Note 

that the assumption that the integrand / is periodic is not really a restriction since 
there are transformations which transform non-periodic functions into periodic ones 
such that the smoothness of the integrand is preserved, see for example 27;.) 

To analyze the properties of a quadrature rule one considers then the worst-case 
^ I error sup^g^^ | J,^ ^,^ f{x) dx — QN,s{f)\, where B-^ denotes some class of functions. 

^ ■ In the classical theory the class e^ of periodic functions has been considered where 

one demands that the absolute values of the Fourier coefficients of the function decay 
sufficiently fast (see [121 [3 [27l \19\). This leads us to the classical measure of the 

quality of lattice rules Pa ~ supjg^s J,^ ^.^ f{x)dx — QN,s{f) , which then for a 
lattice rule with generating vector g ~ (gi, . . . , gg) can also be written as 

h.ezs\{o} 

hg = ( mod N) 

where h — (hi, . . . , hs), h ■ g — ft-igi + • • • + hgQs and \h\ — Wj^i max(l, \h.j\). (Later 
on in this paper we prefer to use the more contemporary notation of reproducing 
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kernel Hilbert spaces, in our case so-called Korobov spaces, (see Section [^75]) . but as 
is well understood (and as is also shown in Section 12.31) the results also apply to the 
classical problem.) 

By averaging over all generating vectors g several existence results for good lattice 
rules which achieve Pa = C(7V""(log7V)"^) have been shown, see p ^ fT3 l iM j [20 l 
[T9l [27] . By a lower bound of Sharygin 26 this convergence is also known to be 
essentially best possible, as he showed that the worst-case error is at least of order 
N^" {log Ny^^. But, except for dimension s = 2, no explicit generating vectors g 
which yield a small worst-case error are known. For s > 3 one relies on computer 
search to find good generating vectors g and many such search algorithms have been 
introduced and analysed, especially recently, see [131 [251 HH IM] ■ 

On the other hand one can of course also use some other quadrature rule QN,sif) = 
X]n=o ^nl{xn) to numerically integrate functions in the class e^. In this case the 
worst-case error in the class e^ for a quadrature rule with weights loq, . . . , clIjv-i and 
points {xq, . . . , xm-i} C [0, 1)*" is given by 

^-1 g27rih-(3;„-3;„,) 

Pa{{xo,...,XN_i})= ^ w„a;„ ^ 1^— . (1.1) 

An explicit construction of such point sets was introduced by Niederreiter, see 
[l6l Theorem 5.3], and is called Kronecker sequence. Here the idea is to choose the 
quadrature points of the form {zk}, fc = 1, 2, . . ., where z is an s-dimensional vector 
of certain irrational numbers (for example one can choose z = (^/pi, . . . , \/Ps) where 
Pi, . . . ,Ps are distinct prime numbers. Depending on the smoothness a certain points 
will be used more than once, see |16j . In practice, problems can occur because of the 
finite precision of computers making it impossible to use points whose coordinates are 
all irrational numbers. 

Another construction of quadrature rules is due to Smolyak |31) and is nowadays 
called sparse grid, see also [8^ . Those quadrature rules are sums over certain products 
of differences of one-dimensional quadrature rules. In principle any one-dimensional 
quadrature rule can be chosen as a basis, leading to different quadrature rules. In 
many cases the weights a;„ of such quadrature rules are not known explicitly but can 
be precomputed. But even if the underlying one-dimensional quadrature rule has only 
positive weights, it is possible that some weights in Smolyak's quadrature rules are 
negative, which can have a negative impact on the stability of the quadrature formula. 
In general, quadrature formulae for which all weights are equal and X]n=o '^" ~ ^^ 
that is, (jjn — N^^ for all n = 0, . . . , iV — 1, are to be preferred. As mentioned above, 
such quadrature rules are called quasi-Monte Carlo rules, to which we now switch for 
the remainder of the paper. 

As the weights for quasi-Monte Carlo rules are given by N~^ the focus lies on 
the choice of the quadrature points. Constructions of quadrature points have been 
introduced with the aim to distribute the points as evenly as possible over the unit 
cube. An explicit construction of well distributed point sets in the unit cube has been 
introduced by Sobol [32]. A similar construction was established by Faure [7] before 
Niederreiter [18] (see also [19]) introduced the general concept of (i, to, s)-nets and 
(t, s)-sequences and the construction scheme of digital (i, to, s)-nets and digital (t, s)- 
sequences. For such point sets it has been shown that the star discrepancy (which is 
a measure of the distribution properties of a point set) is 0{N^^ [log Ny~^) , see [TS] . 
From this result it follows that those point sets yield quasi-Monte Carlo algorithms 
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which achieve a convergence of 0(iV^^(log A^)^*"^) for functions in the class e* for 
aU q; > 2. This resuh holds in the deterministic and randomized setting. 

For smoother functions though, i.e., larger values of a in the class £* , one can ex- 
pect higher order convergence. For example, if the partial derivatives up to order two 
are square integrable then one would expect an integration error of ©(A^^"* (log iV)'^^*)), 
for some c{s) > depending only on s, in the function class e^, and in general, if the 
mixed partial derivatives up to order a/2 exist and are square integrable then one 
would expect an integration error in e^ of ©(^^""(log A^)'^*^'*'"^), for some c{s,a) > 
depending only on s and a. But until now (i, m, s)-nets and (t, s)-sequences have only 
been shown to yield a convergence of at best C'(A^~^(log7V)^''~^) (or 0{N'^^^^) for 
any (5 > if one uses a randomization method called scrambling, see |53]) in e^, even 
if the integrands satisfy stronger smoothness assumptions. 

In this paper we show that a modification of digital (i, m, s)-nets and digital 
{t, s)-sequences introduced by Niederreiter [TSl[Tn] yields point sets which achieve the 
optimal rate of convergence of the worst-case error P2a = C'(iV~^™™'^"''^^(log A^)^*"~^) 
for any integer a > 1 and where d € N is a parameter of the construction which can 
be chosen arbitrarily large. We too use the digital construction scheme introduced 
by Niederreiter [TSl [H] for the construction of {t, m, s)-nets and (t, s)-sequences, but 
our analysis of the worst-case error shows that the t-value does not provide enough 
information about the point set. Hence we generalize the definition of digital {t, m, s)- 
nets and digital (i, s)-sequences to suit our needs. This leads us to the definition of 
digital {t, a, /3, m, s)-nets and digital (i, a, /3, s)-sequences. For a = /? = 1 those def- 
initions reduce to the case introduced by Niederreiter, but are different for a > 1. 
Subsequently we prove that quasi-Monte Carlo rules based on digital (i, a, /3, m, s)- 
nets and digital (i, a, /?, s)-sequences achieve the optimal rate of convergence. Fur- 
ther we give explicit constructions of digital (i,a,min(ck,(i),7Ti, s)-nets and digital 
(i, a, min(a, d), s)-sequences, where d G N is a parameter of the construction which 
can be chosen arbitrarily large. 

Digital {t, 2, 2, m, s)-nets and digital {t, 2, 2, s)-sequences over Zf, (i.e. where a = 
/3 = 2) can also be used for non-periodic function spaces where one uses randomly 
shifted and then folded point sets using the baker's transformation (see [3]). Our 
analysis and error bounds for a = 2 here also apply for the case considered in [3] 
(with different constants though), hence yielding useful constructions also for non- 
periodic function spaces where one uses the baker's transformation. Using a digital 
(i, a, m, s)-net with a scrambling algorithm (see [23]) on the other hand does not 
improve the performance in non-periodic spaces compared to {t, m, s)-nets. 

In the following we summarize some properties of the quadrature rules: 

• The quadrature rules introduced in this paper are equal weight quadrature 
rules which achieve the optimal rate of convergence up to some logA^ fac- 
tors and we show the result for deterministic and randomly digitally shifted 
quadrature rules. The upper bound for the randomized quadrature rules even 
improves upon the best known upper bound (more precisely, the power of the 
log N factor) for lattice rules for the worst-case error in e^ for all dimensions 
s > 2 and even integers a > 2 (compare Corollary 16.51 to Theorem 2 in |20]). 

• The construction of the underlying point set is explicit. 

• They automatically adjust themselves to the optimal rate of convergence in 
the class el^, as long as a is an integer such that a < d, where d is a parameter 
of the construction which can be chosen arbitrarily large. 

• The underlying point set is extensible in the dimension as well as in the 
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number of points, i.e., one can always add some coordinates or points to an 
existing point set such that the quahty of the point set is preserved. 
• TractabiUty and strong tractability resuhs (see [SHj) can be obtained for 
weighted Korobov spaces. 
The outhne of the paper is as follows. In the next section we introduce the nec- 
essary tools, namely Walsh functions, the digital construction scheme upon which 
the construction of the point set is based on and Korobov spaces. Further we also 
introduce the worst-case error in those Korobov spaces and we give a representa- 
tion of this worst-case error for digital nets in terms of the Walsh coefficients of the 
reproducing kernel. In Section [3] we give the definition of digital (t, a, /3, m, s)-nets 
and digital (t, q;,/3, s)-sequences. Further we prove some propagation rules for those 
digital nets and sequences. In Section U] we give explicit constructions of digital 
(t, a,/3,m, s)-nets and digital (i, a, /3, s)-sequences and we prove some upper bounds 
on the t- value. We then show, Section [51 that quasi-Monte Carlo rules based on those 
digital nets and sequences achieve the optimal rate of convergence of the worst-case 
error in the Korobov spaces. The results are based on entirely deterministic point 
sets. Section [6] finally deals with randomly digitally shifted digital (i, a, /3, m, s)-nets 
and (t,a,(3, s)-sequences and we show similar results for the mean square worst-case 
error in the Korobov space for this setting. The Appendix is devoted to the analysis 
of the Walsh coefficients of the Walsh series representation of B2a{\x — y\), where 
-820 is the Bernoulli polynomial of degree 2a. In the last section we give a concrete 
example of a digital (i, a, a, m, s)-net where we compute the i-value by hand. 

2. Preliminaries. In this section we introduce the necessary tools for the anal- 
ysis of the worst-case error and the construction of the point sets. In the following 
let N denote the set of natural numbers and let No denote the set of non-negative 
integers. 

2.1. Walsh functions. In the following we define Walsh functions in base 6 > 2, 
which are the main tool of analyzing the worst-case error. First we give the definition 
for the one-dimensional case. 

Definition 2.1. Let 6 > 2 be an integer and represent fc e No in base b, 

k = Ka^ih"-^ + h Ko with Ki € {0, . . . , 6 - 1}. Further let cjf, = e^'^^/''. Then the 

k-th. Walsh function bwaU : [0, 1) — > {1, w^, . . . ,w|j^^} in base b is given by 

bVfalk[x) = LO^ 

for X S [0, 1) with base b representation x = xib~^ + X2b~^ + ■ ■ ■ (unique in the sense 
that infinitely many of the Xi are different from 6—1). 

Definition 2.2. For dimension s > 2, x — {xi,...,Xs) G [0,1)" and k = 
{ki,...,k,) e Ng we define twal^ : [0,1)'' -^ {1, Wfc, . . . , w^~^} by 

S 

bwalfc(a;) = Y\ bwalfe^. (x^- ) . 

As we will always use Walsh functions in base b we will in the following often write 
wal instead of hwal. 

We introduce some notation. By ® we denote the digit-wise addition modulo b, 
i.e., for X = YlTLw ^i^' ^^^ V = J^'iZw y»^~' ^^ define 

00 
x®y^ ^Zib~% 

i—w 
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where Zi € {0, ... ,b — 1} is given by Zi = Xi + yi (mod b) and let denote the digit- 
wise subtraction modulo b. In the same manner we also define a digit-wise addition 
and digit-wise subtraction for non-negative integers based on the 6-adic expansion. 
For vectors in [0, 1)'* or Ng the operations © and Q are carried out component- wise. 
Throughout the paper we always use base b for the operations © and ©. Further we 
call X € [0, 1) a 6-adic rational if it can be written in a finite base b expansion. 

In the following proposition we summarize some basic properties of Walsh func- 
tions. 

Proposition 2.3. 

1. For all fc, Z G No and all x,y G [0,1), with the restriction that if x,y are not 
b-adic rationals then x <S) y is not allowed to be a b-adic rational, we have 

walfc(a;) • wali{x) = walkQi{x), walfc(x) • walfc(y) = walfc(a; © y). 

2. We have 

1 .1 

walo(x) da; = 1 and / walfe(x) =0 if k > 0. 
la Jo 

3. For all fe,Z G Nq we have the following orthogonality properties: 



I walfc(a;)wali(a;) da; 
^fo.i 



^ 1, ifk = l, 
[0 1)=""^'°^""''""''^""^""' I ^' otherwise. 



4- For any f G £2([0, 1)^) and any cr G [0, 1)* we have 

/ f{xQ)cr)dx^ f{x)dx. 

"'[04)" "'[04)" 

5. For any integer s > 1 the system {wal^ : k = (fci, . . . , fc^), fci, . . . , fc^ > 0} is 
a complete orthonormal system in £2([0, 1)'*)- 
The proofs of 1.-3. are straightforward and for a proof of the remaining items see 
[2] or |33j for more information. 

2.2. The digital construction scheme. The construction of the point set used 
here is based on the digital construction scheme introduced by Niederreiter, see P^ . 

Definition 2.4. Let integers m, s > 1 and & > 2 be given. Let R^ be a 
commutative ring with identity such that |i?{,| = b and let Zj, — {0, . . . , & — 1}. 
Let Ci,...,a e K'""" with Q = (cj,fe,Oi<fcJ<m. Further, let tj^i : Zb ^ Rt for 
I — 0, . . . ,m — 1 and /ij^fe : i?b — > Zf, for j = 1, . . . , s and fc = 1, . . . , to be bijections. 

For n = 0, . . . , 6™ — 1 let n = X^I^o ^li''^)^''^ with all ai{n) G Zf,, be the base 
b digit expansion of n. Let n — {ipo{aQ{n)), . . . ,iljm~i{am~i{n)))'^ and let ijj = 

iVj,!, • ■ • , yj,mV = (^']^ foi" i = 1, • • ■ , s- Then we define Xj^n = iij,iiyj,i)b~^ H h 

l^j,m{,yj,m)b~"^ for j = 1, . . . , s and n = 0, . . . , 6™ — 1 and the n-th point a;„ is then 
given by a;„ = (xi^„, . . . , Xs^n)- The point set {xq, . . . , Xhm_i} is called a digital net 
(over Rb) (with generating matrices Ci, . . . , C^). 

For ?Ti = c» we obtain a sequence {xq, Xi, . . .}, which is called a digital sequence 
(over Rb) (with generating matrices Ci, . . . , C^). 

Niederreiter's concept of a digital (t, m, s)-net and a digital (i, s)-sequence will 
appear as a special case in Section [3l Apart from Section [3] and Section 21 where we 
state the results using Definition 12.41 in the general form, we use only a special case 
of Definiton 12.41 where we assume that 6 is a prime number, we choose Rb the finite 
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field Zf, and the bijections ipi and /ij j. from Z^ to Z^ are all chosen to be the identity 
map. 

We remark that throughout the paper when Walsh functions wal, digit- wise ad- 
dition ©, digit-wise subtraction G or digital nets are used in conjunction with each 
other we always use the same base b for each of those operations. 

2.3. Korobov space. Historically the function class e^ has been used. In this 
paper we use a more contemporary notation by replacing the function class e^ with 
a reproducing kernel Hilbert space Ha called Korobov space. The worst-case error 
expression (|l.ip will almost be the same for both function classes and hence the results 
apply for both cases. 

A reproducing kernel Hilbert space H over [0, 1)'* is a Hilbert space with in- 
ner product (•, •) which allows a function K : [0,1)'' — > M such that K{-,y) E H, 
K[x,y) = K{y,x) and {f,K{-,y)) = f(y) for ah x,y € [0,1)'' and all / G H. For 
more information on reproducing kernel Hilbert spaces see [T], for more information 
on reproducing kernel Hilbert spaces in the context of numerical integration see for 
example [5ll30]. 

The Korobov space "Hq is a reproducing kernel Hilbert space of periodic functions. 
Its reproducing kernel is given by 

where a > 1/2 and \h\ = J|^^-^max(l, \hj\). The inner product in the space Ha is 
given by 



(/,g)„= ^|/ir/»g(/i), (2.1) 



where 



f{h) = / fix)e-'^''^^'^dx 



'[0,1 

1 /? 

are the Fourier coefficients of /. The norm is given by ||/||q = (/, f)a ■ 
Note that for a a natural number and any x G (0, 1) we have 



B2a{x) 



(-l)"+i(2a)! v-^e' 



2-nlhx 



(2.)- ^^ I/.I- 



where B2a is the Bernoulli polynomial of degree 2a. Hence, for a a natural number 
we can write 

J27:±h{xj—yj) \ ^ / fr^ \2a 



2-Ki.h(x-y) (2tt]'^°' 

Kaix, y) = 1 + E -^h^^ = 1 - (-i)"\^^^"(l^ - ^D- (2.2) 
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Then we have 

s 

where x = (cci, . . . ,Xs) and y — (yi, . . . , j/s). Hence the Korobov space is a tensor 
product of one-dimensional reproducing kernel Hilbert spaces. 

Though a > 1/2 can in general be any real number we restrict ourselves to 
integers a > 1 for most of this paper. The bounds on the integration error for Ha 
with a > 1 a real number still apply when one replaces a with \_a\ , as in this case 
the unit ball of Ha given by {/ G Ha ■ \\f\\a < 1} is contained in the unit ball 
{/ e H^a] ■ 11/11 [aj < 1} oiH^al ^^ 11/11 L"J — ll/ll"' Hcncc it follows that integration 
in the space Ha is easier than integration in the space H \a\ ■ 

In general, the worst-case error e{P,H) for multivariate integration in a normed 
space H over [0, 1]" with norm || • |j using a point set P is given by 



e{P,H) = sup 

/e-H.||/||<i 



/ f{x)dx-Qp{f) 

-'[0,1] = 



where Qp{f) = N^^ X^xeP /(^) ^^'^ ^ ^ 1^1 ^^ *^*^ number of points in P. If H is 
a reproducing kernel Hilbert space with reproducing kernel K we will write e{P,K) 
instead of e{P,H). It is known that (see for example [30]) 

e^P, K) = / K{x, y)dxdy--y2 / K{x,,, y) dy + — V K{x^, Xi), 

(2.3) 
where P — {xq, . . . , a;jv_i}. Hence for the Korobov space Ha we obtain 



e2(F,i^„) = -l + ^ ^ if„(a^„,a;^). (2.4) 

nji=0 

Therefore it follows that e^{P,Ka) = P2a and hence our results also apply to the 
classical setting introduced by Korobov [13j. 



It follows from Proposition 12 .31 that Ka can be represented by a Walsh series, i.e., 

let 

Ka{x,y)= ^ rb^a{k,l)walk{x)wali{y), (2.5) 



kjem 



where 



n^aikj)^ / ii:a(a;,y)walfc(a;)wali(y)da;dy. 



'[0,1 

As the kernel Ka is a product of one-dimensional kernels it follows that ri,^a{k,l) 
11^=1 '>^b,a{kj,lj), where k = (fci, . . . ,fcs) and I = (/i, . . . ,Zs) and 



rb,a{k,l)= / i^Q(a:, y)walfe(.T)wal/(y)da;dy. 

Jo Jo 
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For a digital net with generating matrices Ci, . . . , Cg let 2? = 2?(Ci, . . . , Cg) be 
the dual net given by 

P = {fc e Ng \ {0} : Cf fci + ■ ■ • + Cjks = 0}, 

where for k = (fci, . . . , kg) with kj — kj^ + Kjib + ■ ■ ■ we set kj — {kj^, ■ ■ ■ , Kj^m-i)'^ ■ 
Further, for 7^ u C {1, . . . , s} let I?„ = 'D{{Cj)j^u)- We have the following theorem. 
Theorem 2.5. Let Ci,...,Cs e Z™^™ be the generating matrices of a digital 
net Pjjra and let T> denote the dual net. Then for any a > 1/2 the square worst-case 
error in "Hq is given by 

e^{Pb^^,K^)= Y, rtAI^J). 
kdev 

Proof From ((^ and (^3]) it follows that 

e^{Pbm.,Ka)^ -1+ ^ '^b,a{k,l)-r^ ^ walfe (a;)walj (y). 



'52. 



In [5] it was shown that 



L Y- ,.„1...^_/ 1 iffcePuW, 



^ walfc(a;) = I 



^m /_^ '"^ ^ I otherwise. 



Hence we have 



e^ 



fe,;ei>u{o} 



In the following we will show that ri,, a{0,0) = 1 and ri,^a{0,k) = rii^a{k,0) = 
if fc 7^ from which the result then follows. Note that it is enough to show those 
identities for the one dimensional case. We have walo(x) = 1 for all x E [0, 1) and 
hence 



= f walfc(y)d2/+ / J2 1^1"'" / e'^"'^dxe-^^'''ywMy)dy 

= / walk{y)dy. 
Jo 

It now follows from Proposition 12.31 that rf,^Q,(0, 0) — 1 and rii,a{0,k) — for fc > 0. 

The result for ri,^aik, 0) can be obtained in the same manner. Hence the result follows. 

D 

In the following lemma we obtain a formula for the Walsh coefficients rt^a- 
Lemma 2.6. Let b > 2 be an integer and let a > 1/2 be a real number. The 

Walsh coefficients ri,^a{k, I) for fc, / € N are given by 

r. (k D- V ^''■''^^'' 



/iG 
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where Ph,k = /o c ^'^^'"walfe(x) dx. 



Proof. We have 



•^° "^o /iez\{o} 

= I] \h\-^" f e^^'''^walk{x)dx f e-^^^'^^walzly) dy. 
ftez\{o} ° ° 

The result follows. D 

It is difficult to calculate the exact value of rb^k,!) in general, but for our 
purposes it is enough to obtain an upper bound. Note that rf,^a{k, k) is a non-negative 
real number. 

Lemma 2.7. Let b > 2 be an integer and let a > 1/2 be a real number. The 
Walsh coefficients r^^k, I) for k,l G N are bounded by 

|?'6,a(fc,0P < rb,a{k,k)rb,a{l,l)- 



Proof. Using Lemma 12.61 we obtain 

\pKkWhA \ ^ V- \PhA' V- \Ph 



I ti, ni2 / Y^ \PhM\PhM ^ sr^ \i^h,k\ \r^ 

\ftez\{o} ' ' / hez\{o} ' ' hez\{o} ' ' 

= rb,a{k,k)ri,^a{l,l)- 

The result follows. D 

In the following we will write rb^k) instead of rb^a{k, k) and also rb^k) instead 
oirb,a{k,k). 

Lemma 2.8. Let Ci, . . . , C^ G Z™^™ be the generating matrices of a digital net 
Pbrn and let V denote the dual net. Then for any natural number a the worst-case 
error in "Ha is bounded by 



e{Pb^,K^)< Y^^JnAk)- 
kev 

Proof. From Theorem 12.51 and Lemma 12.71 it follows that 

e'{Pb^,K^) < J2 \rbAkJ)\ < [ ^ ^^^^^^(fc^ ) 
k,ieT> \kev ) 

and hence the result follows. D 

For a > 1 a natural number we can write the reproducing kernel in terms of 
Bernoulli polynomials of degree 2a. Then for fc > 1 we have 

^mW - (-1)"+'^S^ / / S2„(|x-y|)walfc(x)waU(y)dxd2/. 



(2a)! h h 
Note that the Bernoulli polynomials of even degree 2a are of the form 

B2o.(x) = co^x^'^ + c„_ia;2("-i) + . . . + cq + cx^''-\ 
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for some rational numbers Cq,, . . . , Cq, c with Cq, c 7^ 0. Let 

1 rl 



hi^) = I I |a;-yPwalfc(x)walfc(y)da;d?/. (2.6) 

Jo JO 

As mentioned above, rb^a{k) is a real number such that rf,^Q,(fc) > for all fc > 1 and 
a > 1/2, hence it follows that for any natural number a we have 

nAk) < -^7^4]- {\Cal2a{k)\ + |Ca-l/2(a-l) (fc) I + ' ' ' + |co/o(fc)| + |c/2a-l|) • 



Using Lemma 18.21 and Lemma 18.51 from the Appendix we obtain the following 
lemma. 

Lemma 2.9. Letb,a eN withb>2. For k e n with k = Kib"'^-'^ ^ ^ Ht.b"'--^ 

where v > I, ki, ■ ■ ■ ,k,i, G {1, . . . ,b — 1} and 1 < a^ < ■ ■ • < ai let qb,a{k) = 
jj-ai a,-„i„(„_„) ^ y/jgy^ Jqj- Q^y natural number a and any natural number b > 2 there 
exists a constant Cb.a > which depends only on b and a such that 

rbAk)<Cl^ql^{k) forallk>l. 



Let now qb,a{^) = 1. For fe = (fci, . . . , ks) € Ng we define qb,a{k) = 0^=1 qb,a{k])- 
We have the following lemma. 

Lemma 2.10. Let m > 1, b > 2 and a > 2 be natural numbers and let X'Jm „ = 
X>„n{l,...,6™- 1}I"I. Then we have 



kev 
< E (1 + ^""'"C'fc,„(a + &-2))^-l"lcl"J^(l + a + 6-2)l"lQL„,n,a(Ci, . ■ • , C,) 

07t«C{l,...,s} 

+ (l + 6-"™a„(a + 6-2)). _i^ 



where Cb.a is the constant from Lemma \2.9\ and where 

^,«,a(Cl>■■•>C's) = Yl IbA^^)- 



keVT 



Proof. Every fe € Nq can be uniquely written in the form k = h + b"^l with 
he {0,..., 6™ - 1}" and I e Ng. Let P^™ = X> n {0, . . . , 6" - 1}^ Then we have 



J2\/'''>Ak)= Y. \jnAb"'l)+ Y E V^^^"(^ + ^'"^)■ 



For the first sum we have 



Y \/rbAb"'l) = -'^+Y y^bAb^'l) = -1 + II \/rbAb"'l) 
By using Lemma 18.81 from the Appendix and Lemma 12.91 we obtain that 



Y yrbAb'^i) = 1 + fe""" Y y^bAi) < 1 + b-'"^Cb,c Y ibAi)- 

1=0 1=1 1 = 1 
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We need to show that J^ili Qb.ai^) < a + b~^. Let / — lib'^^^^ + • • • + l^h^"^^ for some 
V > 1 with 1 < Cy < ■ ■ ■ < ci and li, . . . ,1^ G {1, . . . ,b — 1}. First we consider the 
sum over all those I for which 1 < v < a. This part of the sum is bounded by 

a ooci — 1 c,y_i— 1 a oo 

I/— 1 Cl— iyC2— i^ — 1 Ciy — 1 U — l C—1 

If j^ > a we have qb,a{l) = <lb,a{l') for I — li¥'^~^ + • • • + li,b'^''~^ and where V = I' {I) = 
lih''^~^ + • • ■ + lab'^"~^ . Thus we only need to sum over all /' (i.e. natural numbers 
with exactly a digits) and for given V multiplying it with the number of I which yield 
the same /', which is 6"^°"^ — 1 (and which we bound in the following by 6^°~^). We 
have 

CXD Ci— 1 Cq_i — 1 

(^-ir E E--- E b-^---^-h^--^ 

oo Ci — 1 Ca:-2— 2 

= b-\h-ir Y. E--- E (c.-i-2)6-^---^°- 

ci— a+1 C2— CK Cq — 1— 3 



oo oo 

c=l c=l 

_ 1 



Thus we obtain X^l^i Qb,a{l) < ex + b ^. 
Further we have 

s oo 

E E\AvS^^^= E nE\Av^(v^^> 

heDfcm isNg heV^rri ] = 1 1=0 

where /i = (/ii, . . . ,hs). By using Lemma 18.81 from the Appendix and Lemma 12.91 we 
obtain 

oo CX) 

E ^nJr^) = 1 + &-""a„ Y. 'i^Ai) < i + &""'"a,a(a + &-2). 

Let now < hj < 6™. From Lemma [2]9] we obtain 



\Jn,a{hj + &™0 < Cb^aqb,a{hj + h^l) < Cb,aqb,a{hj)qb,a{l)- 

From above we have X^i^o 1b,a{l) < I + a + b^^ and hence 



^ \/''6,q(^j + ^"0 < '?6,Q(^j)C'b,aE'^^^"'^^) <Cb,a(l+a + 6 '^)qb,a{hj). 



1=0 1=0 
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Thus we obtain 



E Y.y'^bAh + b^l) 



= E E nEv'^^"('^.+^"onEv''M(^"0 

lll^uC{l,...,s}h^eT>',„ ^jeu 1=0 ji^u 1=0 

0#«c{i,...,s} h^^vi„,^^jeu 

where h^ — {hj)j^u- The resuh follows. D 

In [26] it was shown that the square worst-case error for numerical integration 
in the Korobov space can at best be of 0{N~'^°' (log Ny~^), where N is the number 
of quadrature points. Hence Lemma 12.101 shows that it is enough to consider only 
Ql „ „ a(Ci, . . . , Cs) in order to investigate the convergence rate of digitally shifted 
digital nets. 

3. {t, a, (3, m, s)-nets and {t, a, (3, s)-sequences. The t value of a (t, m, s)-net 
is a quality parameter for the distribution properties of the net. A low t value yields 
well distributed point sets and it has been shown, see for example [6l [19] , that a small 
t value also guarantees a small worst-case error for integration in Sobolev spaces for 
which the partial first derivatives are square integrable. 

In the following we will show how the definition of the t value needs to be modified 
in order to obtain faster convergence rates for periodic Sobolev spaces for which the 
partial derivatives up to order a are square integrable. It is the aim of this definition to 
translate the problem of minimizing the worst-case error into an algebraical problem 
concerning the generating matrices. (This definition can therefore also be used in an 
computer search algorithm, where one could for example search for the polynomial 
lattice with the smallest t{a) value which in turn yields a small worst-case error for 
integration of periodic functions.) 

For natural numbers a > 1, Lemma 12.91 suggests to define the following metric 
fit,^a{k,l) = Hb,a{k G I) on Nq which is an extension of the metric introduced in [TT] . 
see also ^5\ (for a = 1 we basically obtain the metric in [171 [25] ). Here Hb.a{0) = 
and for fc G N with k = n^h"-"^^ + • • • + K,ib°'^~^ where 1 < a^ < ■ ■ ■ < ai and 
Ki G {1, . . . , 6 — 1} let ^b,a{k) = ai + - • ■+ai-nin{a,iy)- For a A; G Nq with k — (fci, . . . ,ks) 
let ^b,a{k) = fJ-b.aiki) -f • • • -I- Hb,aiks)- Then we have qb^k) = 5-^^o('^). Hence in 
order to obtain a small worst-case error in the Korobov space "Hq, we need digital nets 
for which min{/ib^c,(A;) : k G P} is large. We can translate this property into a linear 
independence property of the row vectors of the generating matrices Ci, . . . , C^. We 
have the following definition. 

Definition 3.1. Let m, a > 1 be natural numbers, let < /3 < a be a real 
number and let < i < /3to be a natural number. Let Rb be a ring with h elements 
and let Ci, . . . , a G i?^""" with Cj = (c^- 1, . . . , Cj.mV . If for all 1 < «j>^ < • • • < 
ij,i ^ "^j where Q < Vj < m for all j' = 1, . . . , s, with 

«14 H + H,min(i/i,a) H + «s,l H h «s,min(iy, ,q) < Z?™ " * 

the vectors 
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are linearly independent over Ri, then the digital net which has generating matrices 
Ci,...,Cs is called a digital {t,a, (3,m,s)-net over Ri,. Further we call a digital 
(t, a, a, m, s)-net over Rb a digital (i, a, m, s)-net over Ri,. 

If t is the smallest non-negative integer such that the digital net generated by 
Ci, . . . ,Cs is a digital (t, a, /3, ?7i, s)-net, then we call the digital net a strict digital 
(t, a, /3, 771, s)-net or a strict digital (t, a, m, s)-net if a = /3. 

A concrete example of a digital (t, a,/3,m, s)-net, where we also calculate the 
exact i- value by hand, is given in Section [71 

Remark 1. Using duality theory (see [21]) it follows that for every digital 
(i, a, /3, m, s)-net we have minfegp l^b,a(k) > /3m— t and for a strict digital (i, a, 13, m, s)- 
net we have minfegx> l-^b,a{f^) — I^it^ — t+1. Hence digital {t, a, (3, m, s)-nets with high 
quality have a large value of (3m — t. 

Definition 3.2. Let a > 1 and t > be integers and let < /3 < a be a 
real number. Let Rb be a ring with b elements and let Ci,...,Cs G i?|j"^°° with 
0/ ~ (cj,i, Cj_2, • ■ O"^- Further let Cj^m denote the left upper m, x m, submatrix of 
Cj. If for all m > t/f3 the matrices Ci,m, . • . , Cs,m generate a digital (i, a, /?, m, s)- 
net then the digital sequence with generating matrices Ci , . . . , Cs is called a digital 
(t, a, /3, s)-sequence over Rb- Further we call a digital (i, a, a, s)-sequence over Rb a 
digital (t, a, s)-sequence over Rb. 

If t is the smallest non-negative integer such that the digital sequence generated 
by Ci , . . . , Cs is a digital (t, a, j3, s)-sequence, then we call the digital sequence a strict 
digital (t, a, /3, s)-sequence or a strict digital (t, a, s)-sequence if a = /?. 

Remark 2. Note that the definition of a digital (i, l,m, s)-net coincides with 
the definition of a digital (i, to, s)-net and the definition of a digital (t, 1, s)-sequence 
coincides with the definition of a digital (i, s)-sequence as defined by Niederreiter [TS] . 
Further note that the i-value depends on a and /3, i.e., t = i(a, /3) or t = t{a) \ia = (3. 

In the following theorem we establish some propagation rules. 

Theorem 3.3. Let P he a digital {t,a, /3,m,, s)-net over a ring Rb and let S be a 
digital {t, a, /3, s) -sequence over a ring Rb- Then we have: 

(i) P is a digital {t' , a, (3' , m, s)-net for all \ < (3' < j3 and all t < t' < j3'm and 
S is a digital (t' , a, /?', s) -sequence for all 1 < (3' < j3 and all t < t' . 

(ii) P is a digital {t',a' , l3' ,m,s)-net for all 1 < a' < m and S is a digital 
(t' , a' , /3' , s) -sequence for all a' > 1, where (3' = (3imn{a, a')/a and t' = 
[imin(a, a')/a] . 
(Hi) Any digital (t, a, m, s)-net is a digital ([to'/a] , a', to,, s)-net for all I < a' < a 
and every digital [t, a, s)- sequence is a digital {\ta' la\, a' ,s)- sequence for all 
1 < a' < a. 

Proof. Note that it follows from Definition 13.21 that we need to prove the result 
only for digital nets. 

The first part follows trivially. To prove the second part choose an a' such that 
a' > 1. Then choose arbitrary 1 < ij^^j. < • • • < ij^i < m with Q < Vj < m such that 



„min(a,a') 

«l,l^ ^*l,min(i/i,Q')^ ^*s,lH |-*s,min(iy,,Q') < ^(3- 



a 



min(a, a') 



We need to show that the vectors 



are linearly independent over Rb. This is certainly the case as long as 

«1,1 H 1- «l,min(iyi,Q) H 1- is.X ^ V is,min{i^,,a) < /3m - t. 
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Indeed we have 



^1,1 I ■ ■ ■ I ^l,min(i^i,Q) 1 ■ ■ ■ I ^s,l + ■ ■ ■ + ^s,min(z^s,a) 
OL .. . 
< — r-p 7t(*1,1 H 1- «l.min(i/i,Q') H + is,l + ' 

miii(Q;, a') . \ L, J 



a{L's ,a') ) 



< m/3 - 



min(a, a') 



niiii(Q;, a') 



<m(3-t, 



and hence the second part foUows. The third part is just a special case of the second 
part. D 

Remark 3. Note by choosing a' = 1 in part (in) of Theorem 13. 31 it follows that 
digital (t, a, m, s)-nets and digital (i, a, s)-sequences are also well distributed point 
sets if the value of t is small, see [19] . 

4. Explicit constructions of digital {t, a, f3, m, s)-nets and digital {t, a, [3, s)- 

sequences. In this section we show how suitable digital (i, a, (3, m, s)-nets and digital 
{t, a, l3, s)-sequences can be constructed. 

Let d > 1 and let Ci, . . . , C^d be the generating matrices of a digital {t, m, sd)- 
net. Note that many explicit examples of such generating matrices are known, see 
for example [71 [TH HH [31] and the references therein. For the construction of a 
{t, a, /3, m, s)-net any of the above mentioned explicit constructions can be used, but 
as will be shown below the quality of the (t, a, 13, m, s)-net obtained depends on the 
quality of the underlying digital (t, m, sd)-net on which our construction is based on. 

Let Cj — {cjs, ■ ■ ■ , Cj^m)'^ ioT j — 1, . . . , sd, i.e., Cj^i are the row vectors of Cj. 
Now let the matrix Ca be made of the first rows of the matrices C/. 



W :. 



^U-l)d+l,- 



then the second rows of C(j_i)d+i, . . . , Cjd and so on till C^j"' is an m x m matrix, 
, Cj l^y where c, ( = Cu.v with I — {v — j)d + u, 1 < v < m 



i.e.,cf^^ijfl,. 



and {] — \)d < u < jd for I = 1, . . . , ?Ti and j — 1, . . . , s. In the following we 
will show that the matrices C{ , . . . , Ci are the generating matrices of a digital 
(t, a, min(a, d),m, s)-net. 

Theorem 4.1. Let d > 1 be a natural number and let Ci, . . . , C^d be the gen- 
erating matrices of a digital (t' ,m, sd)-net over some ring Ri, with b elements. Let 



(d) 



Ad) 



Cs be defined as above. Then for any a > 1 the matrices C} , . . . , C^ 



(d) 



C- 

are generating matrices of a digital {t, a, min(a, d), m, s)-net over Ri, with 



t = min(a, d) t' + 



s{d — 1) min(a, d) 



Proof. Let C- ' — (cL , . . . , c}- ^)-^ for j = 1, . . . , s and further let the integers 
ii,i, • ■ • , *i,iyi, • • • , *s,i, • • • ,is.vs be such that 1 < ij^i,. < ■ ■ ■ < ij.i < m and 

«ia H ^ ii.min(<.i,a) H ^ isS H h is,min(,y, ,a) < min(a, d)m - t. 

We need to show that the vectors 

Jd) id) (d) (d) 

l,il,l I • ■ • I ^l,ii,„^ 1 ■ • ■ J ^s,is.i 1 ■ • ■ 1 ^s,is,u^ 



Jd) 



M) 



are linearly independent over Rb. For j = 1, . . . ,s let Uj — {cj I ^ , . . . , c^ /. ^}. The 

'(j_i)d+i, ■ ■ ■ , Cjd- For j = 1, . . . , s and 



vectors in the set Uj stem from the matrices C, 
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dj — (j — l)d + I, . . . ,jd let Bd denote the largest index such that (e^ — j)d + dj € 
{ij^v^ ■ ■ ■ 7*j,i} ^iid if for some dj there is no such Cd we set e^ = (basically this 
means e^ is the largest integer such that Cd-,ed- G Uj). 

Let d < a, then we have d{{e(^j-i)d+i — 1)+ + • • • + {ejd — 1)+) + J2i=i ^ — *j,i + 

^ij,min{^j,d) where (x)+ = max(a;,0) and Lj = \{{j -l)d+l < dj < jd : Cd^ > 0}|. 

Hence we have 

Li 

d{{e(j-i)d+i - 1)+ H + {Ejd - 1) + ) + Yl I 

1=1 
= d(e(j_i)d+i H h Cjd) - Ljd + Lj{Lj + l)/2 

did — 1) , , 

><e(j-i)rf+i + --- + ej-rf)- ' ^ ' . (4.1) 

Thus it follows that 

, ^, . , d(d-l) , d(d-l) 
d[ei -\ h Csd) < 2^(^J^ H ^ ^j,min{uj.a)) + s <dm-t + s 

and therefore 

t d-1 , 

ei H \- Csd < m - - + s < m - t . 



Thus it follows from the (t' ,m, sd) -net property of the digital net generated by 
Ci , . . . , Csd that the vectors c[J^ ^, . . . , c[J^ ^ , . . . , c),J^ ^, . . . , 0].^^ ^^ are linearly in- 
dependent. 

Let now d > a. Then we have (i((e(j_i)d+i — 1)+ + • • • + (ejd — 1)+)) + J2i=i ' — 

ijs H h ij\min(i/,,Q) + {d~ a)i J ^^in{u,, a), where again Lj = \{{j - l)d + 1 < dj < 

jd : Cdj > 0}|. Hence we can use inequality (I4.ip again. Note that ii^i-aXa(vi,a) + ■ ■ • + 
*s,min(i/,,a) < rn — tjOL and hence we have 

s 

^(ijaH l-ij,,„i„(^^,Q) + (d-a)ij^„ii„(^^,a)) < am-t + {d-a){m-t/a) = dm~dt/a. 

Thus it follows that 

«(ei H h Csd) < 2^(«i4 H ^ «j,min(z.j,a) + (d " aj«j,min(i/, ,a) ) + S 

J = l 

dt d(d-l) 

<dm h s-!^ 

a 2 

and therefore 

t d-1 , 

ei + • • • + Csd < m h s < m — t . 

a 2 

Thus it follows from the (i', to, s(i)-net property of the digital net generated by 
Ci , . . . , Csd that the vectors c[j^^,...,c{j^^,...,c\l^^,..., c^J^ ^^ are linearly in- 
dependent and hence the result follows. D 

In Section[7]we use this construction method to construct a digital (3, 2, 4, 2)-net 
over 7^2 ■ 
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Note that the construction and Theorem 14. II can easily be extended to {t, a, (3, s)- 
sequences. Indeed, let d > 1 and let Ci, . . . ,Csd be the generating matrices of a 
digital (t, sd)-sequence. Again many explicit generating matrices are known, see for 
example [71 [191 ESI [32]. Let Cj = {cj^i,Cj^2, ■ • 0^ for j = 1, . . . , sd, i.e., Cjj are the 
row vectors of Cj . Now let the matrix C be made of the first rows of the matrices 
C(j_i)£;_|_i, . . . , Cjd, then the second rows of C(j„i)(i+ii • • • : Cjd and so on, i.e., 

j — \^iJ~l)d+l,lT ■ ■ TC.jd.l,C(^j_i)d+1.2i ■ ■ ■ T(^jd,2, ■ ■ ■) ■ 

The following theorem states that the matrices C[ , . . . ,Cs are the generating ma- 
trices of a digital (i, a, min(a, d), s)-sequence. 

Theorem 4.2. Let d > 1 be a natural number and let Ci, . . . , Csd be the gener- 
ating matrices of a digital (t' , sd)- sequence over some ring Rb with b elements. Let 
Ci , . . . ,Cs be defined as above. Then for any a > 1 the matrices C} , . . . ,Cs 
are generating matrices of a digital {t, a, min(a, d), s) -sequence over Rh with 

. / ,x / [ s(d — l)mm( a, d) 
t = min(a, d) t' + ^ ' ^ ' ' 



The last result shows that (i, a, /3, m, s)-nets indeed exist for any < /3 < a 
and for m arbitrarily large. We have even shown that digital {t, a, f3, m, s)-nets exist 
which are extensible in m and s. This can be achieved by using an underlying {t' , sd)- 
sequence which is itself extensible in m and s. If the t' value of the original (t' , m, s)- 
net or (t', s)-sequence is known explicitly then we also know the t value of the digital 
[t, a, j3, m, s)-net or (i, a, /3, s)-sequence. Furthermore it has also been shown how 
such digital nets can be constructed in practise. 

In the following we investigate for which values oft, a, s, b digital (t, a, s)-sequences 
over Zfc exist. We need some further notation (see also [53], Definition 8.2.15). 

Definition 4.3. For given integers s,a> 1 and prime number b let db{s,a) be 
the smallest value of t such that a {t, a, s)-sequence over Zf, exists. We have the 
following bound on db{s,a). 

Corollary 4.4. Let s,a > 1 be integers and b be a prime number. Then we 
have 

(s (6-i)g + b + 
a ( - - 1 - logb 



b 



,, , , ,36-1 (26 + 4)Vs^ „ a(a-l) 

<4(s,a) < c^(g-l)-^3Y~" VP^ - + 2a + s ^ ^ \ 

Proof. The lower bound follows from part (Hi) of Theorem l3.3l bv choosing a' = 1 
and using a lower bound on the t- value for (t, s)-sequences (see |22j). The upper bound 
follows from Theorem 14.21 bv choosing d — a and using Theorem 8.4.4 of [23]. D 

5. A bound on the worst-case error in T-La for digital (i,a,/3,TO, s)-nets 
and digital (t, a, (3, s)-sequences. In this section we prove an upper bound on the 
worst-case error for integration in the Korobov space Ha using digital (i, a, /3, m, s)- 
nets and (i, a, /?, s)-sequences. 

Lemma 5.1. Let a > 2 be a natural number, let b be prime and let Ci, . . . ,Cs € 
jmxm ^g ^^g generating matrices of a digital (t, a, /3, m, s)-net over Zt with m > t/ j3. 
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Then we have 

Qlm,uACu . . . , C,) < 26l-l"6-^"+*(/3™ + 2)l"l"-\ 

where Ql^nu a *■' defined in Lemma \2.1(A 

Proof. We obtain a bound on Qj„ m ^i „, for all other subsets u the bound 
can be obtained using the same arguments. 

We first partition the set I?^„ ,-^ ^i into parts where the highest digits of kj 

are prescribed and we count the number of solutions of Cf A:i + • • • + C^kg = 0. For 
J = 1, . . . , s let now ij^a < ■ • • < ij^i < rn with ij^i > 1. Note that we now allow 
ijj < 1, in which case the contributions of those ij^i are to be ignored. This notation 
is adopted in order to avoid considering many special cases. Now we define 

■^h™,{l,...,s}(*l,l' • ■ • J*l,a; ■ • ■ : *s,li • • ■ j*s,q) 

= {ke 'Dl„.,{i^...,s} ■ kj = [Kj,ib'^-^~^ + ■■■ + Kj,ab'^-"-^ + lj\ with 0<lj < b'^-"-^ 
and 1 < Kj^i < b ior j — I, . . . ,s}, 

where [-J just means that the contributions of ij^i < 1 are to be ignored. Then we 
have 

Qb,m,{l,...,s},Q(Cl, ■ • ■ iCs) 

m *l.a-l— 1 m Js,a-l — 1 17-)* f- ■ ■ ■ w 

E-^ Y^ ^ l-'^fc'",{l,...,s}l'l.l' ■ • ■ /'l.g; ■ ■ ■ > '■s,li ■ ■ ■ 1 '■s,f-^\^ ^ 
2_^ ' ■ ■ Z^ ' ' ' 2^ ^ii.i+---+n.c + ---+is.lH his,c ^ ' 

il,l — 1 ^l,a — 1 is, 1—1 is,a— 1 

Some of the sums above can be empty in which case we just set the corresponding 
summation index ij^i — 0. 

Note that by the {t,a, l3,m, s)-net property we have 

l^6",{l,...,s}(^l,l' • ■ • ' *!,"' • • • ' *s,l' • • ■ ' '^s,a)\ = 

as long as «i^i + • • • + ii,a + • • • + is.i + • • • + is.a < Pm — t. Hence let now < 
ii,i, • ■ • , *s,a < TO be given such that ii^i, . . . ,is,i > 1, ij.a < ■■■ < ij.i < m for 
j = 1, . . . , s and where if ij.i < 1 we set ijj = and ii,i+- • •+ii,Q+- • ■+is,i+- ■ ■+*s,q > 
/3m -— t. We now need to estimate |2?jJ,„ ^ si.(*i,i; ■ • ■ : *i,a7 ■ ■ ■ ^ is,i, • ■ • , *s,a)|, that 
is we need to count the number of fe e I?^,„ ,^ ^1 with kj — [kj^i6'ji^^ + • • • + 

Kj^q6*j.°-i + lj\ such that Cf fci H h C^kg = 0. 

There are at most (6 — I)"'* choices for ki^i, . . . , Ks,q (we write at most because 
if ij^i < 1 then the corresponding Hj^i does not have any effect and therefore need not 
to be included). Let now 1 < ki.i, . . . , Kj^q, < & be given and define 

g = '«14Cvji.i H \- l^l.aC^,ti,^ H ^ '«s,1cIjs,i H ^ Ks,aC^Zs,„ , 

where we set cj", = if / < 1. Further let 

^ ^ ('^141 • • • i'^l,ii.„-li • • • 5 '^s,li • • • i'^s,is,Q-l^- 

Now the task is to count the number of solutions I oi Bl — g. As long as the columns 
of B are linearly independent the number of solutions can at most be 1. By the 
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(t, a, P, m, s)-net property this is certainly the case if (we write {x)^ — niax(a;, 0)) 

< a(ii,Q H his, a) 

< ^m - t, 



that is, as long as 



Let now ii^a 



1l,a 



I ^s.a _ 



Pm — t 



a 



^ s . f> ^ 



pm~t 



Then by considering the rank of the matrix 
B and the dimension of the space of solutions oi Bl = Q it follows the number of 
solutions of Bl ~ g is smaller or equal to 6*1°'' Hs.a,-l{f!m-t)/a\ ^]^^g ^g have 



|2^b'",{l,...,a}(*l,l'---'«l 

r 

(6- 1)"" 



< < 



• ■ ! ^S,l 5 ■ ■ • 7 ^s,a j I 

if E|=i ELi *j,i <(3m-t, 
if Ei=i Er=i «j:i > Pm-t 
andE;=i^..«<^, 
(6 - i)"^5n,„+-+i.,„-L(/9™-t)/"J if J2"j^:^ Er=i «j,! > ;3m - t 

and Ej=i«j,a > ^^• 

We estimate the sum (15.11) now. Let Si be the sum in (|5.ip where ii,i + - • • + Js,q: > 
/3to — i and ii^a + ■ • • + is,a < . For an Z > /Sto — t let v4i(/) denote the number 

of admissible choices of ii,i, . . . , is^a such that I — zi,i + • • • + is, a- Then we have 






Ai(/) 



We have Ai{l) < (^^"1"^) and hence we obtain 



si<{b-ir E 



l=Prn-t+l 



sa — \ J U ~ \ sa — 1 



where the last inequality follows from a result by Matousek [TSJ Lemma 2.18], see also 
[6l Lemma 6]. 

Let 15*2 be the part of (I5.1|l for which ii,i + - • •+«s,a > l3m^t and ii,Q + - • • + «s,q > 
^"'~ , i.e., we have 



rn *l,Q-l— 1 m ^s,n-l— 1 

^..(6-irE-- E •••E-- E WT- 

ill— 1 il,a — 1 ^s.l— 1 is.Q— 1 

Jl,a-2 — 1 m is, 0,-2 — 1 



5-L(/Jm-i)/aJ 



hii,Q-iH his,iH his.Q-l 



< 



m"(6- 1) 



-i)/aj 2^ " 2^ "' ^ "' ^ !,ii,i + ---+ii.„-iH his.l + ---+ii^-J 



^L(/3™-i)/aJ -:1^ -:1^ -:1^ ^^ b^ 

21 1—1 zi a__i — 1 is 1—1 is Q-l — 1 



where in the first line above we have the additional conditions ii.i + ■ • • + is.a > /3m ~t 

From the last inequality and ii,a-i + ■ • • + is,a-i > 



^S.Q ^ 



pm-t 



and ii.Q + • 

ii.„ + ■ ■ ■ + is a for Z = 1, . . . , a — 1 it follows that ii i 



"*l,a-l 



+ «s.l+--- + 
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*s,Q-i > [(/3m — t){l — q:^^)J + 1. Let ^2(0 denote the number of admissible choices 

oiii^i, . . . ,ii^a_i, . . . ,is^i, . . . ,is,a-i such that / = ii^iH hn.Q-iH hi^.iH h 

Zs,Q-i. Note that we have A2{1) < Ct(l-7w7^)- Then we have 

i=[(;3m-t)(l-Q-i)J+l ^ ^ ' 

m''{b-iy°' fj\i0rn-t)/a] ^ [(/3to - t) (1 - a^i )J + s(a - 1) 



^L(/9m-t)/aJ (^l_lj-ly(a-l)ljPm-t+iy s(a - 1) - 1 

where the last inequality follows again from a result by Matousek [15j Lemma 2.18], 
see also [BJ Lemma 6] . Hence we have 

S2 <TO''5""6-/5™+*A(/3w-t)(l-a-i)J +s(a-l) 
~ \ s(a — 1) — 1 

Note that we have Ql^n a fi s}(^i' • ■ • ' ^s) — Si + 82- Let a > 1 and 6 > be 
integers then we have 

Therefore we obtain S*! < 6''"6-^™+*-i(^m-i + 2)""-i and ^2 < 6""6-^™+*m"(^m- 
i + 2)'*("~i)-i. Thus we have 

Qlm,a,{i,...,s}iCi, . . . , a) < 26^"6-^™+*(/3m + 2)^"-\ 

from which the result follows. D 

The following theorem is an immediate consequence of Lemma [2.10l and Lemma f5.ll 
Theorem 5.2. Let b be prime, a > 2 be a natural number and let Ci, . . . ,Cs G 

jinyim ^^ ^^^ generating matrices of a digital {t, a, /3, m, s)-net over Z;, with m > t/ j3. 

Then the worst-case error in the Korobov space "Ha is bounded by 

2 (1 + b-"'"Cb,^(« + b-^) + Cfc.„(l + a + b-^){pm + 2)")^ 
eb,™,a(Ci,...,C3J< 6^™-*(/3m + 2) 

+(l + 6-"™C6,„(a + 6-2)).^_l^ 



where Cb^a > is the constant in Lemma[ 

Remark 4. By the lower bound of Sharygin [26] we have that the worst-case 
error in the Korobov space Ha is at most C'(A^~"(log A^)'*"^). Hence it follows from 
Theorem 15.21 that for a digital {t, a, {3, m, s)-net with (3 > a we must have t — 0{{P — 
a)m). Thus in order to avoid having a i- value which grows with m we added the 
restriction /3 < a in Definition 13. II Further, this also implies that a digital {t, a, /3, s)- 
sequence with t < 00 cannot exist if /3 > a, hence /3 < a is in this case a consequence 
of the definition rather than a restriction. 

Remark 5 . Lemma 12.81 also holds for digital nets which are digitally shifted by 
an arbitrary digital shift cr e [0, 1)'' and hence it follows that Theorem 15.21 also holds 
in a more general form, namely for all digital (i, a, /3, m, s)-net which are digitally 
shifted. 
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Theorem 15.21 shows that we can obtain the optimal convergence rate for natural 
numbers a > 2 by using a digital (t, a, rn, s)-net. The constructions previously pro- 
posed (for example by Sobol, Faure, Niederreiter or Niederreiter-Xing) have only been 
shown to be (t, 1,™, s)-nets and it has been proven that they achieve a convergence 
of the worst-case error of 0{N^^ {log N)^^^). 

We can use Theorem 1 5. 2 1 to obtain the following corollary. 



Corollary 5.3. Let b be prime and let c[ , . . . , Ci € Z^^°° be the generating 
matrices of a digital (i(a), a, min(a, d), s)-segMence S over Zb for any integer a > 1. 
Then for any real a > 1 there is a constant Cj^ ^ ^ > 0, depending only on b, s and a, 
such that the worst-case error in the Korobov space Ha using the first N — b™ points 
of S is bounded by 



JY-mm( [aj ,d) 



e,,^,„(Cr,...,Ci^))<CU&*^LoJ) 



Remark 6. The above corollary shows that digital (i, a,min(a,(i), s)-sequences 
constructed in Section |4] achieve the optimal convergence (apart from maybe some 
logA^ factor) of P2a of 0(A^^^"(log A^)^""^^) as long as a is an integer such that 
1 < a < d. If a > d we obtain a convergence of C'(iV~^'*(log A^)^*"~^). 

6. A bound on the mean square worst-case error in Ha for digital 

{t,a, (3,m,s)-nets and digital (i,a,/3, s)-sequences. To combine the advantages 
of random quadrature points with those of deterministic quadrature points one some- 
times uses a combination of those two methods, see for example [BJ [TUJ [TSl [23]. The 
idea is to use a random element which preserves the essential properties of a deter- 
ministic point set. We call the expectation value of the square worst-case error of 
such randomized point sets the mean square worst-case error. 

6.1. Randomization. In the following we introduce a randomization scheme 
called digital shift (see [SJ [15]). Let Pn = {xq, ■ ■ ■ ,xn-i} ^ [0,1)'* with a;„ = 
{xi^ri, • • ■ , Xs,n) and Xj^n = Xj^n,ib^^ + Xj^n,2b^^ + ■■ ■ ioY n = 0, . . . , N - I and j = 
l,...,s. Let CTj^i, tTj_2, • ■ • <= {0,1} be i.i.d. for j = l,...,s. Then the randomly 
digitally shifted point set Pjv.o- = {^Oi • • • , ^iat-i}, 2;„ ~ (-Zi,„, • • ■ , -Zs,™) using a digital 
shift, is then given by 



<-j,n 



= {^3,n,\ ffi CTj.l)^ ^ + {^3,n,2 © Crj,2)& 



for j = 1, . . . , s and n = 0, . . . , A^ — 1, where a;j>,fc ® crj,„ = Xj^n,k + Cj.n ( mod b) (note 
that all additions of the digits are carried out in the finite field Zb). Subsequently let 
Pat = {xq, . . . ,a;Ar_i} and let Pf^^„ be the digitally shifted point set Pjv using the 
randomization just described. 

6.2. The mean square v^rorst-case error in the Korobov space. In this 
section we will analyze the expectation value of e^(P;v,<T, -f^a), which we denote by 
e'^{PN,Ka) — E,[e'^{Ppf^cr,Ka)], with respect to the random digital shift described 
above. We call e^{PN,Ka) the mean square worst-case error. 

From p.4p and the linearity of the expectation operator we have 



TV-l s 

'^ I I— A r r\ I in. I*-' ifAr r\ ii i — i— 

n,l=aj 



i-{PN,Ka)^E[e^{PN,^,Ka)] = -l + ^ J2 H E[i^„ (z,, „, z,, ,)] ■ 
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In order to compute E[_ftrc(zj „, z^ ;)] we need the following lemma, which, in a 
very similar form, was already shown in [B], Lemma 3. Hence we omit a proof. 

Lemma 6.1. Let xi, a;2 G [0, 1) and let zi,Z2 G [0, 1) he the points obtained after 
applying an i.i.d. random digital shift to xi and X2- Then we have 



Tw\ 1 I \ — r? — Ti J walfe(xi)waU(a;2) if k = I, 
E[wal,(zi)waMz2)] = ^ otherw^se. 



Recall that 



Ka{xi,X2)^ ^ rb,a{k,l)wa\k{xi)wa\i{x2), 
k,i=a 



where 



rb,a{k,l) ^ / / Ka{xi,X2)walk{xi)wali{x2)dxidx2. 
Jo Jo 

Let zi,Z2 be obtained by applying an i.i.d. random digital shift to xi,X2. Using 
Lemma 16.11 and the linearity of expectation we obtain 

oo 
E[Ka(zi,Z2)] = y^ T-fc^g (fc) walfc {Xi )walfc (X2 ) , 
k=0 

where rh,a{k) = rb,a(k,k) and rh,a{^) = 1- 
Therefore we obtain 

-. N-l s oo 

E[e^(PAr^o.,i^a)] = -1 + — 2 ^ J| ^ »-b,Q (fc)walfc (xj>)walfe (a;^-,;) . 

n,l=Q j=l fc=0 

Further we have 

W y^ rfc^a(fc)walfc(a;j,„)walfc(xj,;) = 1 + ^ rh^a(fc)walfe(a;„ Q a;;), 

j=i fc=o fce{o,...,6'"-i}=\{o} 

where we write rh^aik) = rii=i ^fc.al^j) foi' ^ — (fcij---j^s)- We have shown the 
following theorem. 

Theorem 6.2. Let b > 2 be a natural number and let a > 1/2 be a real number. 
Then the mean square worst-case error for integration in the Korobov space Ha using 
the point set P/v randomized by a digital shift is given by 

1 ^"^ 

fcGNgVJO} n,l=0 

In the following we closer investigate the mean square worst-case error for digital 
nets randomized with a digital shift. 

Subsequently we will often write e^ ^ a(C'i, . . . , Cs) to denote the mean square 
worst-case error E[e(Pb">,CT, ^a)], where P;,™ is a digital net with generating matrices 
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Ci, . . . , Cg and 6™ points and Pt'^^a- is the digital net P^™ randomized with a digital 
shift. 

Theorem 6.3. Let m > 1, b be a prime number and a > 1/2 be a real number. 
The mean square worst-case error in the Korobov space %„, using a randomly digitally 
shifted digital net over Zf, with generating matrices Ci, . . . , C^ G Z™'*™ is given by 



p2 
h.m.Oi 



(Ci,...,C,)=^rfc,„(fc). 
kev 

Proof. In [5^ it was shown that 

-^ Y, walfc(a;„ Qxi)^—Yl ^M^n) = < 



1 iffcePujo}, 

otherwise. 



Hence the result follows from Theorem 16.21 D 

Remark 7. Theorem 12.51 and Theorem 16.31 now imply that 



y fcGP Y k,l£V 

i.e. the root mean square worst-case error is always smaller than the worst-case error, 
see also Remark [5l 

Remark [5] and also the above Remark imply that the bounds on the worst-case 
error also hold for the root mean square worst-case error. On the other hand, following 
the proofs for the bound on the worst-case error using the criterium for the root mean 
square worst-case error yields a better bound. We outline the results subsequently. 

Following the proof of Lemma 12.101 we obtain 

kev 0^uC{l,...,s} fcGl5*„^ 



where Cb,a is the constant from Lemma 12.91 and 
Cb,a = Ci 



■\ 



Q-l-l 



(62-fo)-l[](62c_^2(c-l))-l, (61) 



The sum X^fcGC lb ai^) "^^^^ '^^^ ^^ bounded using almost the same arguments 
as in the proof of Lemma f5.ll Doing this one can obtain that for a digital (i, a, /3, to, s)- 
net we have 

J2 IbA^) ^ (2fo)l"l"fo~^^^'"~*)+^(/3TO - t + l)l«l"-i. 

Hence we obtain the following theorem. 

Theorem 6.4. Let b be prime, a > I an integer and let Ci,. . . ,Cs G Z™^™ be 
the generating matrices of a digital (i, a, /3, to,, s)-net over Z;, with to, > t/ j3. Then the 
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mean square worst-case error in the Korohov space "Ha is bounded by 
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;;2 



fc,m,a(Cl, • • ■ tCs) 



< 



1 + 6-2-C2„ + {2bnCl^ + C2 J(/3m - i + 1)" - (1 + 5-^— C^ J 



62(/3m-*)-l(^m-t + l) 



-(1 



f,-2crn^2 



1, 



where Ct, q, > is the constant in Lemma \2.9\ and the constant Ch a > is given by 



We can use Theorem [63] to obtain the fohowing corollary. 

Corollary 6.5. Let b be prime and let C['''\ . . . , Cs"^^ G Z^^°° be the generating 
matrices of a digital (i(a),a, min(a, d), s)-segueKce S over Z(, for any integer a > 1. 
Then for any real a > 1 there is a constant C^ ^ a > 0; depending only on 5, s and a, 
such that the root mean square worst-case error in the Korobov space "Ha using the 
first N — 5™ points of S is bounded by 

eb.m.QlL-i ,---,'^s )^^b,s,aO ATraindal .d) 



66 



JYmin( [a] , 



Remark 8. The above corollary shows that the digital (i, a, min(a, d), s)-sequences 
constructed in Section|4]achieve the optimal convergence of P2a of ©(A^"^" (log Ny°'~-^) 
as long as a is an integer such that 1 < a < d. (This convergence is best possi- 
ble for a = 1 by the lower bound in |26).) If a > d we obtain a convergence of 

Using the construction of Theorem 14.11 or Theorem 14.21 it follows that t{a) also 
depends on the choice of d. Hence choosing a large value of d also increases the 
constant factor 6*(L"J) in Corollarv l5.3l and Corollarv l6.5l 



7. Some examples of digital (i, a, m, s)-nets over Z2. In this section we give 
a simple example to show how the nets described in this paper can be constructed. 
We use the construction method outlined in Section ID 

7.1. Example of a digital (0, 2, m, l)-net over Z2. First we use the so-called 
Hammersley net as the underlying digital net, which is a (0,m, 2)-net over Z2. The 
generating matrices for this net are given by 



Ci = 



/I 


... 0\ 




/O ... 


1\ 


■•• 


■■. 


and C2 — 


•■■ 


• ■■ 


[0 ... 


1) 




V 


... 0) 



(7.1) 



(2) 



Now we use the construction method of Section 2] to construct the matrix C]; 

(2) 
i.e. d = 2 in this case. The first row of C^ is the first row of Ci, the second row of 



-(2) 



(2) 



Ci is the first row of C2 , the third row of C{ is the second row of Ci , the fourth 

(2) 
row of Ci is the second row of C2 and so on. Assume that Ci , C2 are mxm matrices 
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where m is even. Then we obtain 

/I 
... 





c[ 



2) 



1 



1 






... 

yo ... 

So for example if to = 4 we obtain 



c[ 



2) 





1 












1 




/I 








0\ 











1 





1 








Vo 





1 


0/ 




0/ 



(7.2) 



.(2) 



The matrix C\ is of course non-singular and therefore the point set one obtains 
are just equidistant points starting with 0. 

Assume that to, is even. Then the digital net which one obtains from C| is a 
digital (0, 1, m, l)-net over Z2 and, at the same time, it is also a digital (0, 2, to., l)-net. 
Note that using the bound from Theorem 14.11 we obtain a i- value of 1 , but by closer 
investigation using Definition 13.11 one can see that the properties also hold for i = 0. 
Hence the i- value obtained from Theorem 14. II is not necessarily strict even if the value 
of the underlying digital net is strict. 

7.2. Example of a digital (i, 2, 4, 2)-net over Z2. Consider the digital (1, 4, 4)- 
net over Z2 with generating matrices given by Ci , C2 above and 



G.= 



^(2) ■ 



Then C} is given by (17.21) and Cj is given by 



/I 1 1 1\ 

10 1 

11 

\0 1/ 

(2) 



and Ca — 



/O 1 1 0\ 

110 1 

1 

\o 1 oy 



a 



(2) _ 



/I 1 1 1\ 

110 

10 1 

VI 1 ly 



Using the digital construction scheme we obtain the points 

(0 0) (i ^) fi ^) (^ ^) (-^ ^) (^ -^) (^ -^) (ii ^) 



(I 11) f3 l^ ^3 l^ (Z 13>) (_5_ J_) (13 7\ (J_ l^ ( 
\4' 16'' M' S-*' *>8' 4''' Vs' le-*' V16' 16''' V16' S-*' V 16' 2''' V 



16' 



It can be checked that this digital net is a digital (1, 1,4, 2)-net, i.e. a digital 



.(2) 



(2) 



(1,4, 2)-net (the first two rows of Cj; and the first two rows of C2 are linearly 
dependent, so the i- value cannot be when a = 1). 

Now we investigate the t- value when a — 2. First note that Theorem 14.11 yields a 
t- value of 4 for a = 2 (d = s = 2). Further the i- value cannot be 2 in this case: we need 
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to consider all cases where iii +ii.min(iyi,2) +*2,i + *2.min(iy2.2) ^ am — t = 2-4 — 2 = 6 
with < i^i,i^2 l£ 4. But by choosing iii — 12^1 = 2 and ii^2 — *2,2 = 1 we obtain 
the first two rows of C'l and the first two rows of C2 , and as those 4 rows are linearly 
dependent it follows that the i- value cannot be 2. Now let us check whether a i- value 
of 3 is possible: we need to have «i^i+zi^niin(i,j_2)+*2,i+*2, 111111(1^2. 2) — 5, hence i^i, i^2 > 2 
is not possible (because then we would have ii,i + ii,2 + i2,i + ^2,2 > 2 + 1 + 2 + 1 > 5). 

(2) 
Further the conditions are satisfied if either z^i = or j/2 = as the matrices C^ and 

(2) 
6*2 are non-singular. If i^i > 2 then ii^i > 3 and ii.2 > 2 and hence ii^i + 11^2 ^ 5 

and we can only get ii^i + ii_inin(j/i,2) + *2.i + *2,min(i'2,2) < 5 if z^2 = 0. Hence if either 

i/i > 2 or z/2 > 2 the properties are also satisfied. Thus we are left with the following 

three cases: (i^i, 1^2) = (1, 1), (j^i, ^^2) = (1, 2) and (z^i, 1^2) = (2, 1). 

(2) (2) 

Now let 1^1 = i/2 = 1 . Then we need to take one row of each matrix C^ and Cj 
such that the sum of their row indices is smaller or equal to 5 and check whether those 
two rows are linearly independent. It can be checked that this is always the case: let 
Cj — {cji,cj2,cj^^,cj,^), i.e. Cj^k denotes the fc-th row of CJ . Then the pairs of 
vectors (ci.fc,C2,i) where fc + / < 5 are always linearly independent for all admissible 
choices of k and / (i.e. ci^k 7^ C2.1). 

(2) 
Consider now i^i = 1 and 1^2 = 2, i.e. we take one row from C^ and two rows 
(2) 
from 6*2 such that the sum of the row indices does not exceed 5. Note that 12,2 has 

to be 1 otherwise 12,1 + 12,2 > 5 and ii^i cannot even be 1. As ii^i > 1 and ^2,1 > 2 

the only choices left are ii^i — 1 and 12,1 =2,3 and ii^i = 2 and «2,i = 2. So we need 

to check whether the triplets (ci^i, C2,i, 02,2), (ci,i, C2,i, 02,3) and (ci.2, C2,i, C2.2) are all 

linearly independent, which upon inspection can be seen to be the case. 

The case v = 2 and z/ = 1 can also be checked as the previous case. In this case 
all the relevant sets of vectors are also always linearly independent, hence a i- value of 
3 is possible for a = 2, i.e. the digital net above is a (strict) digital (3, 2, 4, 2)-net. 

The classical t- value (i.e. a = 1) of this digital net is not as good as for example 
the t- value of the Hammersley net (which is 0). On the other hand it can be checked 
that for a = 2 the t- value of the Hammersley net where to = 4 is 4 and hence for this 
case it is worse than the i- value of the digital net constructed above. 

As a last example let us consider the Hammersley net again for arbitrary m > 1, 
i.e. with the m x m generating matrices given by (j7.ll) . As for example the first 
row of Ci and the last row of C2 are the same (and therefore linearly dependent) we 
must have /3m — t < m + 1 for all a > 1 (for a = 1 we can still choose /3 — 1 and 
t = and hence the Hammersley net achieves the optimal i-value, but for a > 1 we 
have seen in Section |4] that there are better constructions). It is sensible to choose /? 
such that we can have a t- value which is independent of to, (for example this is the 
case when one considers sequences and which is also the motivation for introducing 
those parameters; for digital nets it would of course also make sense to just state the 
value of l3m ~ t and to. instead of i,/3 and m). This means that /? < 1, and as /? 
indicates the convergence rate one can obtain it follows that one cannot expect to 
obtain a convergence rate beyond (6'")"^+'' (for an arbitrary small 6 > 0) when using 
a Hammersley net. 

8. Appendix: Some lemmas. We need the following lemmas. 

Lemma 8.1. Let j>l,a>0,b>2 and < w, u < 6" with u ^ v. Then we have 

f.(u+l)/b'' /•(u+l)/b 



I / \x-yyAx<ly 



/6" A/6" ' " ' 6«0-+2)(j + l)Q- + 2) 
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and 

(«+i)/6" /.(t-+i)/6" 2j! ^^-' |u-i;p-2' 






(i- 201(2/ + 2)!- 



Proof. We have 



/ \x — y\-^dxdy— / / |a; — yP dxdy 

m/&» Ju/b'^ Ja Jo 

1 <•! 



^^LL ^^-y^'^'^^'y- 



We divide the last double integral in two parts, we have 

1 |.i f^ f^ f^ f^ 

/ \x — y\-' dxdy — / {y — xY dxdy + / {x — yY dx dy. 

Jo Jo Jo Jo Jy 

We calculate the first part and obtain 



/ / iv - xy <^x dy = ——- I y-'+My 
Jo Jo J + 1 Jo 



(.? + !)(.? + 2) 
and the second part is given by 

I I {x- yy dxdy = -^ I {I- yy+' dy ^ ,.^^.^..^^. 
Jo Jy .? + 1 Jo U + l)U + 2) 

Hence we have 

nl ,.1 2 

\x — vP dxdw = -. 



|.(«+l)/6" |.(t;+l)/6" ^.l/b" |.(|«-t;|+l)/&" 

/ / |a; — yj-' dxdy = / / |a; — y^ dxdy 

Ju/b<^ Jv/b" Jo Jlu-vl/b" 

1 pi /-lu-vl + l 



For the second part we have 

/ Ix-yPdxdy = / / 

lu/b" Jv/b-^ Jo J\u-v\/b 

1 /•! /-{u-vl + l 

'^ , 

' -' \u—v 

where now |m — u| > 1. We have 

nl i-lu-vl + l 1 r-1 

I / (x-yydxdy=--- {{\u-v\ + l-yy+'-{\u-v\-yy+')dy 

J\u-v[ .? + -L Jo 

_ 2\u - v\^+^ - i\u -v\ + iy+^ - {\u -v\- iy+^ 

" (J + !)(.?+ 2) ■ 

The result follows by simplifying the sum in the numerator. D 

Lemma 8.2. Let k > 1 be given by k = Kai-ib"-'^^^ + • • • + Ka^-ib"-"^^ for some 
J^ > I7 /^ai-i, • • • J i^a^-i G {1, ...,&— 1} and 1 < a,y < • • • < fli. For any even 
< .7 < 2i^ we have Ij{k) = 0. 
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Proof. The result for J = follows from Proposition 12.31 and (|2.6p . It was shown 
in [S], Appendix A, that 



oo 6—1 



c— 1 r— 1 ^ ' 



and hence 

(oc b-l . \ -^ 

c— 1 r— 1 ^ ' / 



°° ' ^ wal^fcc,^i(7/) - wal^fcc,-i(a;) 



CXJ , J 0—1 

ci 

Let 

Then we have 

v4fe(ci, . . . ,c-,) 



JjCiH hCj LLA^ Q-2-Ki.T/b _ -^ 

Ci,...,C, — 1 i — lT — 1 



[^ /^ -pT>r^ waUc,-i(y) - waUc.-i(x) — -— 



ci,...,c,=l 



We have 



-py>-^ wal^bc,-i(y) - wal^j,c,-i(a:) 

iiZ^ Q-2iriT/b _ I 

i=l T=l 

Ti,...,Tj = li=l uC{l....,j} i£u i^u 

= E nCe-'^'^^/"-!)-^ E (-l)'"'^alc„,.(2/)walc,,„.„,j,„,.(x), 

Ti,...,Tj=l i = l uC{l,...,j} 

where Cu.t = X]ie« ''"«^^' '^ ^'^"^ hence 

Afc(ci,...,Cj) 

= E t[(^-''''^'^'-^r" E (-1)'"' 

ri,...,rj=l i=l mC{1,....j} 

/•I /•! 



walc„,^ (y)walc{i j.,\„_^ (x)walfe(a;)waU (y) da; dy 

,...,rj— 1 i— 1 

E (-1)'"' / walc„,,efc(2/)d?; / wa^i_...,,^j\„_,efc(a;) dx. 



/o ./o 

6-1 J 



nC{l,...j} 
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Note that ii ly > j /2 we either have C^^t ® fc 7^ or C{i j}\„ ,- Q fc 7^ and 
hence Ak{ci, . . . , Cj) ~ 0. The resuh now foUows. D 

Let ap{n) — X]ft=i h^- It is known that 



ap{n) = Y^ 



" Bh pi 



^^^Mip + l-hy. 



p+l-h 



where Bq, Bi, . . . are the Bernouhi numbers (in particular, Bq = I, Bi — —1/2 and 
B2 = 1/6). 

Lemma 8.3. Let b>2, l<d<a, k — n^^ib'^^^ + • • • + kq where k^-i € 
{1, . . . , 6 — 1}, Kd-2i ■ ■ ■ T 1^0 G {0, . . . , 6 — 1}, m — raa-ib""^^ + • • • + Jtiq and n = 
Ua-ib"^^^ + • • • + TiQ. Then we have 

b^-a fc"-i /, ^ \ b" 



J2 E ^Minem)/bn=b'''-'(l 



2 Q27T±K.d-l/b _l J 2 

n— m—n-{-l 



b^-a &"-! /-, , ^ fe« 



J2 J2 ("^-«)walfc((nem)/6")==63"-2d ji 



n— r?i— n+1 

onrf 

t"-2 6"-l 



6 2sin (K(i-i7r/6)/ 6 



E E (™-n) = J(6'"-fo") 



n— m— n+1 

Proof. In order to obtain a formula for the first sum, let m' = uia-ib""^^ + • • • + 
ma-d+ib''-''+\ m" = rua-d-ib''-'^-^ + • • • + mo, n' = Ua-ib''-^ + ■■■ + na-d+i&""''+^ 
and n" = na-d-ib'^~'^~^ + ■ ■ ■+no. First consider the case where m' > n' and arbitrary 
m",n". We have 

b-l 6-1 

\^ \^ l^27rl(Ko(na-l-"ia-l)H hKd-l(na-d-ma-d))/'' _ Q 

«a-d=0 m„_d=0 

as X]m=o e^'^^'*™/'' = for all K = 1, . . . , 6— 1. Thus we only need to consider the case 
where m' ^ n' , for which case we have 

2iri(Ko(na-l-ma-l)H hKd-i (na-d-nia-d))/fc _ „27riKd-l(nia-d-na-d)/fc 

This part is now given by 

lyd-1 y^ y^ y^ y^ ^27riKd-l(na-d-nia-d)/fc /g 2) 

n" —0 m"—0 ^a-d— r?ia-d— 

where we have the additional assumption rua-db"''^ + rn" > Ua-db"'''^ + n" . First 



ja-d_^ja-d_^ 6-1 6-1 
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consider the case where rUa-d > ^a-d- This part of (|8.2p is given by 

6""''-16""''-l &-2 6-1 

^d-l y^ y^ y^ y^ g27riKd_i(n„_d-m„_d)/& 

n"— m"—0 rin — d— rria-d— ^a-d + 1 

6-2 6-1 

_ Ld-lj2(a-d) \ ^ \ ^ 2TTiKd-i{na-d~ma,-d)/b 

"a-ci=0 ma-d=na-d + l 

u2a—d 

g2-7riKd-l/f' _ 1 ■ 

Now consider the case where iria-d = na-d- In this case we have the assumption that 
m" > n" and hence this part of (|8.2p is given by 



^' E E 1 - ^ (fe^"-'^ - ^>i 



n"— m"— 7i" + l 

Thus ([5^ is given by 

1.2a — d 1 



ri2'7riKc;_i/?J 



T+ic^ 



2a-d 



b") 



and the first resuh follows. 

For the second sum let again m' = 77ia-i6°~"'^ + ••• + nia-d+ib'^'^'^^^ , rn" = 
TOa-d-i&°"''"^ + • • • + Too, "' = na-ife""^ + • • • + Ua-d+ib''^'^^'^ and also n" = 
na-d-iV''^^^ + • • • + riQ. First consider the case where m' > n' and arbitrary m", n" 
We have 

6-1 6-1 

y^ y^ ("^j _ J2)e^'^^^''°^"''~^~™""^-'"' |-td-l("a-d-'na-d))/b 

no-d=0 m„_d=0 
6-1 6-1 



= E E ('"a-d - n„_d)e2-i( 



l(Ko(na-l— mo-l)H |-Kd-l("a-d— ma-d))/6 

""a — a '"a — aj^ 

na-d=Onia-d=0 

0, 



as E™io e^'^i'^™/'' = for ah K = 1, . . . , fe - 1. 

Thus we are left with the case where m' = n' . We have 

g27ri(Ko(n„_i-ma_i)H hKd-l ("a-d-mo~d))/6 _ g27riKd_ i (m„_d-na-d)/b 



Hence this part is given by 



6-1 6-1 6""''-16""''-l 

^'''' E E E E (m" -n" + b''-'\ma-d-na-d))c^''''''-''-"-''-"'-''^/\ 

na-d—^'''ii'a-d—0 n"— 7n''—0 

(8.3) 
where we have the additional assumption nia-db"' "^ + to" > Ua-db"" "^ + n" . First 
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consider the case where rUa-d > ^a-d- This part of (|8.3p is given by 



Q<7ia-d<^na-d<b m" ,7l"—0 

b~2 b-1 

Ij3a-2rf 

2sin^(K(j_i7r/6) 

Now consider the case where rria-d = na-d- In this case we have the assumption that 
m" > n" and hence this part of (|8.3p is given by 



ja-d_2 b'-''-! 



f'" Y Y {m"-n")^-{b'^--''^-b^-''). 



n"— rn"— n" + l 



6 



(This resuh can be obtained using (|8.ip . see the proof of the third part below.) Thus 
is given by 

1.5a — 2d ud 



2sin (K(;_i7r/6) 6 

and the second resuh fohows. 

The third resuh can easily be verified by using (jS.ip . Indeed we have 



b°-2 6"-l 



b°-2 6°-l-n 6"-2 



6°-2 



^ ^(m-n)^^ ^ ^^^^^(6«_„)^l^((5«-„)2_(5«-r^)). 



n— m— n+1 



n— m— 1 



The last sum can be written as i X]n=i("-^ — n) = ^(o'2(&° + 1) — cri(6° + 1)) and by 
using (|8.ip again the result follows. D 

Lemma 8.4. Let j > Q, v > I, I < a^ < ■ ■ ■ < ai < a, k ^ Kai-i&°^^^ + • ■ • + 
Ka,^-ib'^''~^ where Kqi-i, • ■ • , i^a„-i G {l, ...,&— 1}. T/ien we have 



fc"-2 b"-l 



n— m— n+1 



a + l)(j+2) 



and 



b'^-2 b°-l 



y^ y^ (m — n)^walfc((n to)/&°) 



n— 771— n+1 



< Cfe ,60'+2)°-2(ai + ---+o„i„(„,rj/21)) 



for some constant Cb.j > which is independent of v, a and ai, . . . ,a^. 
Proof. We have 



b"-2 b"-l 



b"-l 



^ ^ (m - ny -$](&"- n)n^ = 6V,(fe'^) - <7,+^{b^), 



Tl— 771 — 71 + 1 
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and by using (IS.ip it follows that 
6V,(6'^)-a,+ir) 

°(j + l)(j+2)' 

from which the first part follows as i3o = 1- 

For J = 0, 1 the second part immediately follows from Lemma 18.31 Let now j > 2 
and assume the result holds for all j — 1, . . . , 1, 0. 

Let m — nia-ib"^^^ + • • • + m-o and n — Ua-ib""^^ + • • • + uq. In order to obtain a 
bound on 



y^ y^ (^ _ jj')Jg27ri(K„j_i(no_„j-m<,_„jH |-k„„ _i (n„_„,, -m„_„„ ))/6 

n— ?n— n+1 



(8.4) 



we first sum over the digits ma-ai and Ua-at ■ 

Let m' = ma-ib^-^ + - ■ •+TOa-ai+i&"""^+\ n' = n„_i&''-i+- • •+n,_,,+i&«-«i+i, 
m" = ma-ai-ib°'~"'^^^ + • • • + mo and n" — Ua-m-ib"'^'^^^^ + • • • + no- We consider 
two cases, namely where m' > n' and where m' — n' . 

For m' — n' we either have rria-ai > "-a-ai or rria-ai = na-ai and ?7i" > n", as 
m > n. First let ma_aj > Ua-ai- We have fe"^!^^ choices for m' = n' and the sum 
over the digits ma-ai, "-a-ai with iria-ai > "-a-ai can be written as one sum so that 
the part of (|8.4p where m' = n' is given by 



^ai-l 






n''— m"— r— 1 
bo-^l-lb"-"!-! b-1 

^^"'"^ H H ^(6-T)(T6"-"i+m"-n")-'' 

n"—0 va"—fd T—1 

for some constant C^' ■ > which only depends on b and j. Hence this part satisfies 
the bound. Now let ma~ai = JT-a-ai, then we have m" > n" and hence the part of 
(|8.4p where m' = n' and iria-ai = "-a-ai is given by 

4^ -^ J + 1 J + 2 

„"=0 m"=ri" + l ^-^ '^-^ ' 

where the inequality was already obtained in the first part of this proof. Hence also 
this part satisfies the bound. 
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Now we consider the part of (|8.4p where m' > n' . We have 



6-1 



^ (m' - n' + b^-^^ima-a, - Ua-a,) + m" - „")Je2'^^''"i-i(""-"i-""-"i)/'' 



6-1 



b{m' -n' + to" - n")^' + ^(6 - r)[e-2'^''^"i-i^/''(m' - n' + rfe"-"! + to" - n")^' 
6(to' - n' + to" - n")^' + ^ r j (to' - n' + to" - n")^-"5"(°-"i)£„, 



(8.5) 



M = 



where 



6-1 



Eu =^(6-r)[e^27rlK„,_ir/fc^u^g2,ri«„,_ir/fc(^_^)«j^ 

T = l 

It can be checked that Eq = —b and i?i = 0. Hence (|8.5p is given by 
^ P") (to' - n' + to" - n")^-"5"(°-°i)K, 



«=2 



and hence the resuh foUows from the induction assumption or the first part. D 

Lemma 8.5. Let k > 1 be given by k = Kai-ib"-'^^^ + • • • + '«a„-i^"'' "'^ for some 
v > 1, 1 < flix < • • • < ai and Ka^-i, ■ ■ ■ , Ka„-i G {1, ■ ■ ■ ,b — 1}. Then for j > 1 we 
have 

for some constant Cbj > which depends only on b and j . 

Proof. Let k = Ka-ib""^^ + • • • + kq, where now a = ai, u — Ua-ib'^^^ + •••+1*0 
and V = Va-ib"'^^ + ■ ■ ■ + vq. Tlien we have 

^sO^)^ / / |a;-2/Pwalfc(x)walfc(y)da;dy 
Jo Jo 



6°-16°-l ^(u+l)/6" <.(d+1)/6" 

«=0 «=0 "'"/''° "'"/''° 



For u = u we have e'i^^iM-^<^-i~^,^-i)+-+'^c.-i{uo-v„)) ^ ^ ^gj^^g Lemma 1511 it follows 
that this part in the above sum is given by 



6-0+i)(j + l)(j+2)- 
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Hence it remains to calculate 

^ ^ g2^i(Ko(«„-i-t-„-i)+-+K„_i(uo-^o)) / / \x-yydxdy 

M=0 "=o •lu/b'^ Jv/b'^ 

h"-2 6°-l ., U/2J I |,_2-i 

-^ ^ 6'^0'+2) Z^ (7--2i)!(2i + 2)! 



4j! 



i— u— V— ii+1 



2^ r-,- - 9iV('9i4-9M 2^ 2^ 



5aO+2) ,^ („■ _2i)!(2i + 2)! -^ ^-^ (v-uY^-'J 



where we used Lemma 18.11 The absolute value of the inner double sum can now be 
bounded using Lemma [8.41 and hence the result follows. D 

Lemma 8.6. Let b > 2 be an integer and let a > 1/2 be a real number. Then we 
have 



^rfc,„(A:) = 2C(2a), 
fe=i 

w/iereC(2a) = Er=i'^"^"- 

Proof. Let h e Z\ {0} and let fh{x) = e^^^"". The Walsh coefficients A(fc) 
of the function fh arc then given by fh{k) = L fh{x)'walk{x) dx. It follows that 

\fh{k)\'^ — \Ph.k\'^, where /3/is. was defined in Lemma [2.61 Using Parseval's equality 
we obtain 

oo oo ^1 ^1 

E l/3''.'»l' = E 1^(^)1' = / \hix)\'dx^ ldx = L 

1,-1 ,,_i "'0 Jo 



fc=l fc=l 

Hence we have 

T.^Uk)= E u^El/^'.^l'- E U^-2C(2a). 
fc=i ftez\{o} ' ' fe=i /iez\{o} ' ' 

The result follows. D 

Lemma 8.7. Let k ~ Kai-ib°-^~^ + • • • + Ka^-i^"""^ with 1 < a^ < ■ ■ ■ < ai ana 
let Kai~i, . . . , Hat,-! G {1, ■ . ■ ,b — 1} . Then 

P?,,K„j_ib<'i-i+...+K„„_ib"^-i - 2^ (27ri)'' l-l- 






{2^±Y 1^1 /i. 



Proof. First we consider /c = Ka-i&°^^ with Ka-i G {l,---:^ ^ !}• Let x = 
^ + fr + ■ ■ ■ , then we have walfc(x) = e'^-^^'^o.-ixa/b^ Note that walfc(a;) is constant in 
the intervals [u/6", {u + l)/b°-) for < w < 6". Let li == Ma_ife°"^ H h wq- Then for 
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any h E Z\ {0} we have 



\ ^ ^2TTiKa-iUa/b / 




^-2^ihx ^^ 




u/b" 




6"-l 


-27Tih 


(«+l)/&° _ ^-2-Klhu/b'' 


u=0 




— 2'!T±h 


. 9^ih/(," b-1 




6-1 



y . . . y e2-i'=-i 



'"o/b„-27ri/i(«„_i/bH I-mo/6") 



27r/ii 

110— ''ia-l— 

1 _ „-27ri/i/6" ^-1 *>-! , &-! 

_ ^ "^ V^ 2Triuo{K.a-l/b-h/b'') V^ -27ri«i/i/6""^ V^ -2irltt„_i/i/6 

2Trh± ^ ^ " ' Z^ 

ttO— Ul— tta-l— 

Let now /i € ^\{0} and let /i — hc-ib'^~^ + - ■ ■ + ho and set he = hc+i = • • • = 0. li h > 
we assume that hi S {0, . . . , fe— 1} and if /i < we assume that /li G {—b+ 1, . . . , 0} 
for all i > 0. li ho ^ then E«^!i=oC"^'''"°"'''^'' == and hence /3,j^«^_^,,a-i = 0. 
li ha ^ then E!1^!i=o e^^"'""""'''/'' = &• In general, if for an < i < a - 1 we 
have hi ^ Q then Ph.Ka-iht^-^ — 0. Further, if /ij = for < i < a — 1 then we 
also have /3h.K„_ib»-i ~ 0. Hence, in order to obtain jSh.Ka^ib'^-'^ 7^ we must have 
/iQ = • • • = ha-2 = and Ka-i ^ ^a-i = (mod 6). In this case we have 

-| p— 27ri/ia — 1/6 

Ph.Ka-lb'^-'^ — r> ■ ■ ^ > 

■ ZTTAll 

where /i == ha-ib°'~^ + /ia&° + • • • with ft-a-i = ^a-i (mod b). We can also write 

6(1 - c-2'^iVb) 



/^/ifc^-l^Ka-lb"-! 



27ri/i 



with ft, e Z such that ft = Ka-i (mod 6). 

We can interpret jS^i^ — L e^'^^^'^^walk(x) dx as the Fourier coefficients of the 
k-th Walsh function, hence it follows that 

walfc(a;) = ^/3,,,fce2-i'*^ 

Let now k = Ka^-ib"^^^^ + ■ ■ ■ + Ha^-ib"'"^^ for some 1 < a^ < • • • < ai. Then we 
have 

wal«,^^_^{,ai-i+...+^^^_^fca,.-i(a;) 
= wal^^^_^fcai-i(x) • • • wal^^^ _^f,»--i(a;) 

/iiGZ /ii,GZ 

/il,...,/ii,eZ 

On the other hand we have 

/iez\{o} 



Explicit constructions of quasi-Monte Carlo rules for periodic integrands 35 

On comparing the last two equations we obtain that /3^^^ _i6"i-ih vk^ _i6«>'-i — 

if either 6°i~i |^ or ft. ^ Kai-i (modfo^i^^). Now let /i £ Z such that b°-^^'^\h and 
h = Koi-i (mod 6°!^^). Then we have 

hl,...,h„eS,h;Sjt„j_l ( mod 6) 



h=hifa"l~-^H ^-h^h'^v- 



^ (27ri)^li 



5^ JLl - Q'^T^ihi/b 



/il ,. .. ,/ijy GZ,h^ =/^Q _-L ( mod b) 



i2niy l\ hi 



and the result follows. D 

Lemma 8.8. For k > 1, b >2, m > I and a > 1/2 we have 



Proof. First note that /3?i_„^ _jb™+«i-i^ ^^^ _j(,T,i+a„-i = if fe™ /fft. Further it 

follows from the previous lemma that 

and hence by Lemma 12.61 we have 

r, (hh'"''\ - \^ |/3feb"',fcb'"| _ ,-2am V^ \PhM _ ,-2am (u\ 

hez\{o} ' ' /iGZ\{o} ' ' 

The result follows . D 
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